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Abstract — We characterize stability regions of two-user fading 
Gaussian multiple access (MAC) and broadcast (BC) networks 
with centralized scheduling. The data to be transmitted to the 
users is encoded into codewords of fixed length. The rates of the 
codewords used are restricted to a fixed set of finite cardinality. 
With successive decoding and interference cancellation at the 
receivers, we find the set of arrival rates that can be stabilized 
over the MAC and BC networks. In MAC and BC networks with 
average power constraints, we observe that the duality property 
that relates the MAC and BC information theoretic capacity 
regions extend to their stability regions as well. In MAC and 
BC networks with peak power constraints, the union of stability 
regions of dual MAC networks is found to be strictly contained 
in the BC stability region. 

I. Introduction 

Traditionally, information theorists and the networking com- 
munity have taken two different approaches to the study 
of communication networks. The information theoretic ap- 
proach explores the capacity region [1] of a network by 
optimally designing the encoding and decoding schemes while 
abstracting the higher layers to infinitely backlogged queues. 
The networking approach on the other hand explores the 
stability region [2] of a network by optimizing the scheduling 
and routing algorithms while abstracting the physical layer 
encoding and decoding schemes to server processes. 

Both networking and information theory perspectives have 
lead to many important results. For example, the networking 
approach has found stability-optimal scheduling and routing 
algorithms [3], [4], [5] such as the maximum weight match 
scheduling and dynamic back-pressure routing algorithms. 
One of the remarkable successes of information theory, on 
the other hand, has been the successful characterization of the 
capacity regions of the Gaussian multiple access (MAC) and 
broadcast (BC) channels. Information theory has also revealed 
an elegant duality relationship between the MAC and BC 
capacity regions [6]. 

Both networking and information theory perspectives have 
their own limitations. A common limitation of the network- 
ing research has been that sophisticated encoding/decoding 
techniques introduced by information theory, e.g. successive 
decoding, are overlooked. The utility of information theoretic 
results on the other hand is limited by the numerous underlying 
ideal assumptions such as infinite backlog, continuous rate 
adaptation and infinitely long codewords. In practice com- 
munication systems are designed with a finite set of fixed 



rate codebooks. The rate adaptation algorithm picks one of 
the available rates depending on channel and interference 
conditions, for example [7] [8]. 

In this work, we combine the strengths of information 
theory and networking approaches to study wireless broadcast 
and multiple access networks. Motivated by the networking 
approach, we study the stability regions of wireless multiple 
access and broadcast networks when only limited (finite) sets 
of rates are available to the physical layer Motivated by 
the information theoretic approach, we allow the physical 
layer to use some of the optimal multiuser encoding and 
decoding schemes for the MAC and BC. In particular we 
allow independent encoding and successive decoding with 
successive interference cancellation for both these networks. 
Note that independent encoding and successive decoding suf- 
fice to achieve all points within the information theoretic 
capacity region of the Gaussian MAC and BC. Also, motivated 
by the information theoretic duality of the MAC and BC 
capacity regions we explore the existence of similar duality 
relationships for the stability regions of the multi-rate MAC 
and BC networks. 

There has been extensive work on the design of scheduling 
and routing algorithms aware of both the channel state and 
queue state [9], [10]. [3] proves the modified largest-weighted 
delay first (MLWDF) scheduling policy for broadcast networks 
without interference cancellation at the receivers using the 
concept of Lyapunov drift. [4] extends this to the optimal 
dynamic back-pressure scheduling and routing algorithm to 
generic multi-hop networks with peak power constraints. [5] 
extends these scheduling policies to networks with average 
power constraints and attempts to find distributed implemen- 
tations of these policies in networks. The dynamic back 
pressure algorithm is adapted to cooperative networks in [11]. 
Associated with the problem of finding optimal scheduling 
policies is the problem of determining network capacity and 
stability regions of communication networks. The stability 
and network capacity regions of general networks with peak 
power constraints have been found (Refer [12]). [13] and 
[14] characterize stability regions of multiple access channels 
with a probabilistic multi-packet reception model (MPR) with 
distributed scheduling policies - the ALOHA protocol in 
particular [15] studies capacity and stability regions of certain 
large regular wireless networks with the MPR model. Early 
work in combining information theory with the networking 



approach can be found in [16], [17]. Much like our work, [18] 
incorporates queues into the information theoretic framework 
and explores the equivalence of stability regions of a multiple 
access networks and their ergodic capacity regions; however, 
unlike our work, it does not assume a finite number of 
codebooks. 

The rest of the paper is organized as follows. Section 
defines and discusses stability and networks capacity regions 
of networks. Section presents the system model and the 
problem is defined in Section|lV] Stability regions of MAC and 
BC networks with peak power constraints are found in Section 
[Vl Section [VTl describes techniques on finding stability regions 
of MAC and BC networks with average power constraints. 
Section IVllI presents duality results between MAC and BC 
stability regions. The stability characterization of these multi- 
rate networks leads to the problem of design of optimal 
codebooks discussed in Section IVlllI Section |IX] concludes 
with some interesting directions of future work. 

II. Background 
We review two equivalent definitions of the stability region. 

A. Stability region of a network 

A stationary stochastic process Q{t) is defined to be stable 

if 

lim limsupPr[Q(i) > x] = 0. 

A communication network is defined to be stable if the 
stochastic processes representing its individual queue states 
are stable. The stability region of a network is defined to be 
the closure of the set of all arrival rate matrices that can be 
stabilized by some centralized scheduling and routing policy 
with complete knowledge of the arrival statistics, the queue 
state and other parameters of the network such as channel 
conditions, distance between nodes etc. 

B. Equivalent definition based on infinitely backlogged nodes 

It is shown in [15] that under very mild assumptions that are 
true for most communication networks (i.e. interference does 
not increase capacity), the stability region is equivalent to the 
network capacity region (distinct from information theoretic 
capacity region) defined as follows. 

The network layer capacity region is defined as the set of 
all average departure rate matrices that can be achieved in 
a network assuming that all nodes are infinitely backlogged 
taking into consideration all possible centralized scheduling 
routing and resource allocation strategies in the network with 
possibly complete knowledge of the parameters of the network 
such as channel conditions, distance between nodes etc. 

It must be noted that the information theoretic capacity 
region is different from the network capacity region; the 
definition of the former includes information theoretic as- 
sumptions such as asymptotically long codewords and un- 
countable number of codebooks to choose from, whereas the 
latter definition is more general. The equivalence of stability 
region and network capacity region of a network established 



in [15] is important because the latter is formulated under 
the assumption of infinitely backlogged nodes. The infinite 
backlog model eliminates the dependence on instantaneous 
queue states from the stability optimal scheduUng and routing 
algorithms. 

III. The system model 

We shall study the stability regions of multi-rate Gaussian 
multiple access (MAC) and broadcast networks (BC) (Figure 
[TJ. We first describe the generic system model that is common 
to both the models of multiple access and broadcast networks 
in this work. For ease of exposition, we deal with 2 user 
MAC and BC networks in this work. However the techniques 
described and the duality results extend to general N-user 
MAC and BC networks. 

The channel is a Rayleigh fading channel with additive 
white Gaussian noise at the receiver/s (one receiver in the 
MAC network and multiple receivers in a BC network). The 
noise power at all receivers in all the networks considered are 
assumed to be unity. The channel fading is represented using 
the gain vector h{t) = {hi{t), h2{t)), each component of the 
vector corresponding to a particular link in the network. Note 
that since stability regions and network capacity regions are 
to be identical to each other, we use the network capacity 
model of infinite backlog for our analysis. The following are 
important aspects of the system model. 

A. Encoding 

The transmitters transmit a single fixed-length codeword 
per time slot (i.e. the length of one time slot is normalized 
to the length of the codeword). The codeword is chosen 
from a fixed finite set of Gaussian codebooks. Let sets 
7^^ = {0,R}),R\...R\j_^},i = 1,2 indicate the rates of 
the codebooks available for transmission corresponding to 
user i (rate indicating that the user does not transmit). 
Therefore, the transmitted rate vector in a particular time slot 
lies in TZi x TZ2- For example, in a MAC network with both 
users having a fixed codebook of rate Rq, the rate vector 
that can be transmitted in a particular time slot belongs to 
{(O,0),(i?o,0)(0,i?o)(i?o,i?o)} (indicated by x symbol in 
Figure |2]i. 

B. Slow Fading - Long coherence times 

The channel h{t) is ergodic and remains constant during a 
time slot. The coherence time of the channel is assumed to be 
long enough (and therefore a time slot is long enough) so that 
codewords of sufficient length can be transmitted in a given 
channel state so as to guarantee information theoretic relia- 
bility. The slow fading assumption is reasonable for wireless 
systems with low mobility. 

C. Decoding 

The decoding strategy at the receivers is restricted to succes- 
sive decoding and interference cancellation from the received 
signal. The receiver can decode one complete codeword at a 
time and cancel the corresponding interference and then pro- 
ceed to decode the next codeword. This successive decoding 



procedure is carried on until the receiver decodes his own 
codeword. 

For example, in a two user multiple access network, there 
are only two different ways of decoding corresponding to two 
possible decoding orders. The codeword corresponding to user 

1 can be decoded first, treating the complete codeword corre- 
sponding to user 2 as noise. This is followed by cancellation 
of interference of user 1 and user 2's codeword is subsequently 
decoded. The receiver could alternately choose to decode using 
the reversed order i.e user 2 can be decoded first, followed 
by interference cancellation and then user 1 can be decoded. 
Consider a fixed-codebook MAC network i.e a MAC network 
of two users with each user having a fixed codebook of rate 
Rq. As indicated earlier, one of the four points indicated by 
' X ' symbol in Figure |2] are transmitted. Consider a channel 
state h and a power allocation {Pi{t), P2{t)) such that the 
information theoretic capacity region of the MAC network at 
this channel state is represented by pentagon 1 in Figure |2] 
Then, a codeword of rate (ri,r2) that can be successfully 
transmitted in this channel state and power allocation is 
achievable with one of two possible decoding orders. If rate 
pair (ri,r2) is achievable with user 1 decoded first followed 
by user 2, it must belong to the rectangle with vertices A, B, 
E, O. Similarly, if rate pair (ri, can be achieved with user 

2 decoded first, then it must belong to the rectangle F, C, D, O. 
Therefore, in time slot t, the rate pairs that can be transmitted 
reliably over the network with this power allocation lie in the 
'L-shaped' region (traced in thick lines in Fig. |2]i. It should be 
noted that in Figure|2] although (Rq, Rq) lies in the pentagon 1 
- the pentagon representing the MAC capacity region, it cannot 
be achieved with either decoding order and is therefore not 
achievable at a time slot t with this particular power allocation. 

Note that in a BC network, in a particular time slot, 
all points in the BC capacity region can be achieved with 
successive interference and decoding. The optimal decoding 
orders at both receivers depends on the channel state. 

IV. Problem Definition 

We aim to study stability regions of multiple access and 
broadcast networks. The discussion in section |ll] implies that 
we can assume that each node is infinitely backlogged and 
study the set of achievable departure rates. We categorize the 
problems we solve into the following: 

1) To find the set of supported rate pairs (ri, available 
for a scheduler at a given channel state. 

2) Use the results of part 1 to characterize the stability 
regions of multiple access and broadcast networks. 

3) Explore duality relationships between stability regions 
of multiple access and broadcast networks. 

The rate of transmitted codewords are limited by constraints 
on the power of the transmitted codewords. Power constraints 
usually appear in the following forms: 

« Fixed Power constraints - These are power constraints 
of the the form P{t) ~ P,V<. Analysis of systems with 
this form of power constraints is useful, particularly when 
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Fig. 1. 2-user Multiple Access and Broadcast Networks 



Pentagon 1 : supports {(0, 0), (i?o, 0), (0, Rq)} 
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Pentagon 2 
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Fig. 2. Capacity region and set of supported rates of a MAC at particular 
channel state (i.e hi = h2 = 1) for two different power allocations viz. 
(10, 10) and (10,4) 



the systems are being studied for a short period of time. 
However, this model may not be realistic for analysis of 
certain systems over a long duration. 

• Peak Power constraints - Power constraints of the form 
P{t) ^ P,yt are called peak power constraints and 
model limits on the instantaneous power of the trans- 
mitted signal imposed by the limits of the power source 
and the amplifiers in the system. 

• Average Power constraints - Average power constraints 
are of the form E[P{t)] < P and are used to ensure a 
certain minimum battery life for the users. 

Note that a fixed power constraint is a stronger constraint 
than a peak power constraint, which is in turn stronger than an 
average power constraint. Also note that in our system model, 
a fixed power constraint is not equivalent to a peak power 
constraint since better rates can be achieved by transmitting 
less than the maximum possible power For example, in Figure 
|2] pentagon 2 represents the MAC capacity region with a 
power allocation strictly smaller than that which achieves 
pentagon 1 at the same channel state. We have already argued 
in Section UlI-CI that {Rq,Ro) cannot be transmitted with the 
power allocation corresponding to pentagon 1. Since {Rq, Rq) 
lies in the 'L-shaped' region associated with pentagon 2 
(plotted with broken lines in Figure |2]i it is achievable with 
the corresponding power allocation. Therefore reducing the 
power allocation from that which achieves pentagon 1 to that 
corresponding to pentagon 2 makes {Rq, Rq) achievable in the 
same channel state. In this work, we focus on MAC and BC 
networks with peak and average power constraints. 



In the remaining part of this work, we use the following 
notation: 

F{h,-)^\og{l + h'-) 
X — (xii X2) = (/ii, ^12)^ ~ power of channel fade 



X 



Hadamard reciprocal of x 



V. Stability Region of MAC and BC Networks 
WITH Peak Power constraints 

The stability region of a general single hop network with a 
finite number of channel states and peak power constraint is 
known ([19]) to be: 



r-^pr(5)Co(C(5)) 

s 



(1) 



where pr(S') denotes the probability with which the channel 
takes state S, C{S) denotes the set of all rate vectors that 
can be transmitted in the network reliably, when the channel 
state is S and Co{X) denotes the convex huU of set X. For 
the Gaussian MAC and BC networks, while the Rayleigh 
fading channel vector h{t) — {hi{t), h2{t)) takes values 
over a continuum, a finite state space is naturally defined 
by the partitioning of the channel space based on the rate 
pairs that can be supported. Let C ~ 2^1^^^ denote the set 
of all subsets of TZi x TZ2 with cardinality K = \C\. Let 
Ci G C,i — 1,2- --K denote all possible sets of supported 
rates (or equivalently all possible elements of 2'^'-^''^^). Then 
define the state space S"; , i = 1 , 2 . . . i^T to be 

S, = {{hi, h2) : e the complete set of 

supported rate vectors, if the channel is {hi, ^,2)} 

Note that the set of supported rates is unique to each channel 
state and equation [T] can now be applied to obtain stability 
regions. It must be noted that Si as defined above may be 
empty for certain i. For example. Si = (fiM Ci = {{Rm, Rm)} 
because whenever reliable transmission is possible with rate 
vector {Rm, Rm), reliable transmission is also possible at rates 
(0, Rm), {Rm, 0) and (0, 0). Therefore {{Rm, Rm)} can never 
be the complete set of supported rates for any channel state 
(/ii, /12) (in both the MAC and the BC networks). 

The definition of Si can also be understood using Figure 
|2] Given (/ii,/i2) and powers Pi < Pi,P2 < P2, an L- 
shaped region similar to Figure |2] can be drawn. Let us 
represent by A{hi, /i2), the union of all such L-shaped regions 
over {Pi,P2) pairs satisfying the power constraints. Then 
(/ii,/i2) £ S'i if and only if the set of all rate vectors 
from TZi x 7^2 that lie in A{hi,h2) is precisely Ci i.e 

n A(/ii,/l2) = C,;. 

A. Characterization of Si for a MAC network 

A rate vector (ri,r2) can be achieved in a MAC network 
with one of the two decoding orders. If the rate vector is 
achieved with user 1 decoded before user 2, then conditions 
I2] and |3] below have to be satisfied. 
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) > ri (3) 



where P2min is chosen such that F{h2, ^^j'" ) = ^2- Con- 
dition [3] indicates that the codeword corresponding to user 1 
can be decoded when user 2 is transmitting at the minimum 
possible power(so as to cause minimum interference). Note 
that this condition can be rewritten as 

F{hi,^)>ri 

Similarly, for achieving (ri,r2) with user 2 being decoded 
first, the channel state {hi,h2) must satisfy F{hi,^) > ri 
and F{h2, ^)>r2. Define ^^1^2(5") as follows 

Hi^2{S) — {{hi,h2) : all elements of S can be achieved 

using decoding order (1,2) } 

where S G 2^^^^^. Using conditions |2] and [3] it can be shown 
that 



Hi^2{S)^{{hi,h2):hi> 



2 ^ (2'-i - 1)2-= 



Pi 



2^2 _ 1 



P2 



V(ri,r2) e 5} 



F{h2,^)^\0g{l + h'^P2)>r2 



(2) 



Similarly, 

H2^i{S) = {{hi,h2) : all elements of S can be achieved 

using decoding order (2,1) } 

can be characterized. We can now express Si as follows. 

= i?2^i(c,) u i/i^2(c,) n {Hi^2{C!) u H2^i{c^)r 

where A'^ represents the complement of set A. The above 
equation is merely another form of the definition of Si 
suggesting that it is the region of the channel state space 
where all rate vectors in C, are supported and the rate vectors 
outside d are not supported. Note that the characterization of 
Hi^2{S) indicates that the non-empty channel states Si can 
be represented as rectangular regions in the (/ii,/i2) plane 
(see, for example. Figure O The stability region of the MAC 
network can then be calculated using equation [T] 

Example : Fixed codebook MAC network 
Consider a MAC network with a single fixed codebook of rate 
i?o at both users, i.e TZ = TZi = TZ2 — {0,Ro}. Therefore 
we have 2** = 16 possible subsets of TZ x TZ. However, 
it is easy to see that, except for five of these subsets, the 
channel partitions corresponding to all other subsets are empty. 
These five subsets, represented as Ci, i — 1,2, 3, 4, 5 and their 
corresponding channel partitions Si are described below. 

1) 5i : This represents the state where the channel gains 
are so low that transmission is not possible on either 
channel, i.e i^(/i,,P,) < Ro,i = 1,2 with Ci = {(0,0)} 

2) 52 : Transmission is possible only on channel 1 i.e 
C2 = {(i?o,0),(0,0)} and therefore F{hi,Pi) > Rq, 

F{h2,P2) < Ro 
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Fig. 3. Optimal mapping of channel state (h\,h2) to rate sets d in the 
MAC network with peak power constraints 



3) 53 : Transmission on channel 2 — {{Q,Rq), (0,0)} 
and F(hi,Pi) < Ro, F(/i2, A) > Ro 

4) 54 : Transmission is possible individually on both chan- 
nels AND simultaneously on both channels i.e C4 = 
7^ X 7^ = {(0, 0), {Ro, 0), (0, i?o)(i?o, i?o)}. For (0, Ro) 
and (i?o,0), we require 

F{h^,P^) >Ra,i=l,2 

{Rq,Ro) can be achieved with two possible decoding 
orders. If user 1 is decoded first, user 2 transmits at the 
minimum power required to decode a codeword of rate 
Rq so that the user causes the minimum possible inter- 
ference. Therefore, the condition for successful decoding 
of {Rq, Rq) with user 1 decoded first can be written as 

Pi 



Similarly, if user 2 is decoded first, the condition for 
successful decoding is. F{h2, j^^^a^^ — — ) > -Ro where 
F{hi,Pi,n,n) = -Ro and F{h2, P2min) = Ro Therefore, 
the condition for h to lie in 6*4 can be described by the 
following equations 

F{h,,P,)> Rn,i = 1,2 

F{hi, ) > Ro or F{h2, ^) > Ro 
5) ^5 : Transmission is possible individually in both chan- 
nels, but simultaneous transmission is not possible, i.e 

C5 -{(O,0),(i?o,0),(0,i?o)} and 

Fih,,P,) >Ro,i^l,2 
Pi 



F{hi, 

F{h2, 



2Ro) 

P2 
2Ro] 



< Rq 



< Ro 



Figure [3] shows the partitioning Si of the {hi, /12) space and 
the sets of supported rates Ci corresponding to those partitions. 



The stability region of this network can be found ana- 
lytically from the definition in ([T]i and is plotted in Fig. |6] 
for Rq — 1 and power constraints Pi and P2 satisfying 
Pi + P2 = 2. (The variance of the channel fade is assumed 
to be unity). It can be shown that as long as the probabilities 
pr(S'i) are greater than zero, the stability region is a pentagon. 

B. Characterization of Si for a BC network 

In a BC network, the channel state imposes a decoding order 
i.e if hi > /i2, then both users decode the codeword corre- 
sponding to user 2 first; node 1 then decodes the codeword 
corresponding to user 1 after canceling the interference from 
the codeword corresponding to user 2. Therefore, for hi > /i2, 
the condition for achieving codeword (ri , r2) can be expressed 
as 

F{hi,P)>ri 



F{h2, 



P-Plr. 



-) > r2 



1 + h2Plr 

where is chosen such that F{h2, Pimin) = i"!- Define 

Hi^2{S) and H2.i{S) as follows 

Hi,2{S) — {{hi, /12) : hi > /i2, all rate vectors of set S can 
be achieved in the broadcast channel } 

H2,i{S) — {{hi, /12) : hi < /i2, all rate vectors of set S can 
be achieved in the broadcast channel } 

Hi 2{S) can be characterized as 

^^1,2(5") = {{hi,h2) : hi > h2,Ar,,r,X-^ < (P, P)^ , V(ri , Ts) 
where 

A 

^ri,r-2 - ^ (•2'~i - 1)2''2 (2'^^ - 1) 

and as defined earlier, x is the power of the channel fade. 
H2.i{S) can be similarly characterized. Therefore, 

5, - i?2,l(C,) U Hi^2{Q) n {i?l,2(Q) U -ffl,2(Cf)} 

Note that the characterization of iJi 2(5') indicates that the 
boundaries between different 5*^ in the {hi, ^,2) plane can be 
expressed as linear equations in (for example see Figure 
S. 

Example : Fixed codebook broadcast network 

Consider a two user broadcast network with fixed codebooks 

of rate Rq and the transmit power governed by P{t) < P. 

Similar to the MAC case, the channel state is divided into 5 

non-empty bins listed below. 

In the description below Pimin and P2min take values such 



that log(l + hlPirmn) = -Ro and log(l + hlP2.. 



Ro 



1) Si : Ci = {{0,0)}, ^ F{h„P) < Ro,i = 1,2 

2) ^2 : C2 = {(0,0),(0,Po)}, F{hi,P) < Ro and 
F{h2,P) > Ro 

3) S3 : C3 = {(0,0),(Po,0)}, ^ F{h2,P) < Ro and 
F{hi,P) > Ro 
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Fig. 4. Optimal mapping of channel state {hi,h2) to rate sets d in the 
BC network with peak power constraints 

4) ^4 : = {(O,0)(i?o,0)(0,i?o)(i?o,i?o)} 

F(/i„P) > i?o,i= 1,2 

^(^i'TTW^)>^o if /ii</i2 

and 

F{h,,^—^^)>Ro if /ii>/i2 

the latter two equations implying that the weaker user 
is able to decode the codeword, with the stronger user 
transmitting at the minimum possible power, so as to 
cause minimum interference. 

5) 55: C5 = {(0,0)(i?o,0)(0,i?o)} 

F{hr,P) >Ro,i^ 1,2, < Ro and 

Figure |4] shows the partitioning Si of the (/ii, /12) space and 
corresponding rate sets Ci in the BC network. Figure |6] has a 
plot of the stability region of a broadcast network (which is 
also a pentagon) with P = 2 and i?o = 1- 

VI. Stability regions of MAC and BC Networks 
WITH Average Power Constraints 

The problem in finding stability regions of networks with 
average power constraints is a little more involved because 
of the following reason. In networks with average power 
constraints, at any given instant of time, all the available rates 
can be transmitted by pumping sufficient power. Therefore, 
the problem also involves finding the optimal power allocation 
for each channel state, and therefore the simple expression in 
equation [T] cannot be applied. We therefore use techniques of 
convex optimization to find the stability region ([20] uses the 
technique in a similar context). 

The stability region is equivalent to the set of all average 
departure rates assuming infinite backlog. A scheduler can 
therefore use time-division to stabilize any convex combina- 
tion of rates in the stability region implying that the stability 
region is convex. Since the boundary of a convex region 



can completely be described by its tangents, the boundary 
of the stability region can be found by finding all tangents 
of the region. This is equivalent to maximizing < WjR > 
,V0 < w < 1, under the given power constraint, where 
R — {Ri , R2) represents the average rate vector in the stability 
region, w — {w,l — w) represents the slope of the tangent, 
and < x,y > represents the dot product of vectors x and y. 
Note that this formulation is due to the convex nature of the 
stability region and is independent of the form of the power 
constraint. In general, the optimal rate allocation strategy may 
involve transmitting different rate vectors (possibly randomly 
with different probabilities) at a particular channel state. For 
example, for fixed codebook MAC networks with peak power 
constraints, it was optimal to multiplex between {Rq,0) and 
(0, Ro) when the channel state was in ^5. However in net- 
works with average power constraints, if the channel comes 
from an infinite state space, the optimal strategy transmits a 
single unique rate vector at a given channel state. We state 
this formally below. 

Lemma 6.1: In MAC and BC networks with average power 
constraints, if the components of h take values from a con- 
tinuous state space, the optimal rate allocation strategy that 
maximizes < w,R > associates with a channel state h a 
unique rate vector R{h) G TZi x TZ2- And therefore, the 
optimization problem may be expressed as 

^ max < w, R{h) > f{h)dh 

R{h)e'JZixTZ2 Jh 

s.t / Pi(h)f{h)dh < Pi,i e {set of transmitting nodes} 

Jh 

and R{h) can be reliably transmitted at channel state h with 
powers Pi{h) 

Proof: Proof is presented in Appendix U ■ 
The problem reduces to one of finding the optimal R{h) for 
a given w, so that the power constraint is satisfied. We now 
treat the MAC and BC networks separately. 

A. MAC Network 

Consider a MAC network with the power constraints at the 
two nodes described by E[Pt{t)] < Pi,i = 1,2. In the MAC 
network, for a given channel state, rate vector can be achieved 
with different transmit powers based on the decoding order 
Let us denote two possible decoding orders of a two-user 
MAC network by by 11 = {711,712}- tti represents the order 
of decoding user 1 first and then user 2, and 712, vice-versa. 
Let Pj ' {r,h),i = 1,2 denote the transmit powers of the users 
required to transmit rate vector r at channel state h when the 
decoding order is fixed at tTj . The solution to the optimal rate 
and power allocation problem is given by 

(P*(x),7r*(x)) = arg max {<w,r> 
^>iiP^if,x)-^2P^ir,x)} 



with the Lagrangian muhipUers ki and K2 chosen so that the 
power constraints of the system are met. It can be shown (using 
equations similar to |2] and [3] in IV-Al i that 

2r2(2''i - 1) 



hi 



{T'^ - 1) 



hi 



where r = (r'i,r2) 



(2''i - 1) 



and 



^^^^ 2'-i(2-^-l) 



hi 



where r — {ri,r2) The set of rate vectors (ri,r2) G 
^ divides the channel state space x into different regions 
5" - . = {X : R*{x) = (^1*, '^2), ^*(X)}- It can be shown that 



(l'l.I'2) 



- ?'1,''2 


— > and 




flir*)- 
flir*)- 


flir) f^iri) 
f?ir) fHr^*) 


-flir) ) 
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fm flir*) 
fUr) fiir*) 


-fKrl \ 
-flir) J 



with flir) = 2-^(2'-i - 1), /i(r) = (2'^^ - 1), /^(r) = (2^1 - 
1) and = 2''i(2''2 - 1) 

Note that the boundaries between various rate regions are a 
set of linear equations in 

The maximum value of < w, r > can the be found as 



< w,f >— 



<w,fl> prih e S^-. ) 



For example, in a fixed codebook system with Rq being 
the rate of the codebook, the optimal scheduler would choose, 
when the channel state is h, the maximum of the following : 

1) (No transmission) 

2) wRo — Ki- — ri — - (Transmission only along channel 1) 

3) (1 — w)Ro — K2- — — - (Transmission along channel 2) 



4) 



(2«" - 1) (2«" - 1)(2^°) 
Ra - f^i r? '^2- 



hj 



hi 



5) 



( (i?o, Ro) is transmitted with user 2 decoded first ) 

(2««-l) (2^"-l)(2^o) 

i?0 — K2 — Ki- 



hi 



hi 



((i?o,i?o) is transmitted with user 1 decoded first) 
Figure |5] shows the mapping of the optimal rate to the 
channel state for w — 0.5. (Note: for w — 0.5, ki — K2 
because of symmetry). Note that the boundaries between these 




Fig. 5. Mapping of cliannel state to optimal rate vector in MAC network witli 
average power constraints witli ui = 0.5, (Lagrangian multipliers assumed to 



be Ki = K2 = A;) with 70 : 



Rq 



and 71 = 2^0 - 1 



the Hadamard reciprocal of 



channel states are linear in x ^ 
X- 

Note that any increase in Ki has the effect of reducing the 
corresponding average power consumed E[Pi (t)] . In fact, Ki = 
corresponds to no power constraint and Ki ^ 00 has the 
effect of E[Piit)] 0. This can be observed in Figure|5]where 
the effect of increasing k will be to increase the size of the 
region transmitting (0, 0) suggesting that the power consumed 
reduces. This monotonic behavior is exploited in simulations 
in locating (k^;. Kj) - the Lagrangian multipliers which ensure 
that the power constraints are satisfied. (Note that in the power 
constraint the inequality can be replaced by equality). 

The stability region of a fixed codebook MAC networks 
with Rq = 1 and average power constraints Pi and P2 
satisfying Pi + P2 = 2 are found by performing the above 
optimization for all < w < land presented in Figure |71 

B. BC network 

In the broadcast network, the channel state imposes a 
decoding order Therefore, for a given channel state x> the 
optimal rate allocation is given by 

R*ix) = E^rg max {< w.f> — KP(r, /i)} 

with K chosen so that the transmitter power constraint is met. 
Since, for a given channel state, the decoding order is fixed, 
it can be shown that 



Pir,h) 



2''i-l I (2''2-l)2'' 



if hi > /i2 
if h2 > hi 



The rate vectors fl e TZi x TZ2,i = 1,2...L divide the 
channel state space into regions 

={X: irl,r*2)=R*ix)} 

. We can write 



where 



and 



S(rt,r',~> = {X : XI > X2, {rl,r*^) = i?*(x)} 



= {X : XI < X2, (rt,r*) = R*ix)} 
. It can be shown that, for k = 1,2 

where w =< w.f* — r> and 

aI^) ^ (2'-i - 2''S (2''2* - l)2'^i - (2'^^ - l)2''i) 

The optimal rate allocation for a broadcast network with 
a fixed codebook of rate i?o chooses, at channel state h, the 
maximum of the following; 



1) 



(2-'^'' - 1) (2^'' - 1) 
max{0, wi?o — ^ ri ■, (1 ^ w;)_Ro ~ ri } 



2) if /ii < /i2 



, (2^"-l) ^ (2^°-l)(2^°) , 



/l2 



and if h\ > /12, the following quantity 

.(2^«-l) (2^0 - 1) (2^0). 



hl 



hi 



The optimal rate allocation partitions the channel state in a 
manner similar to the MAC ( Figure |5]l, for hi = h2 and 
w = 0.5. 

The stability region of a fixed codebook BC network with 
i?o = 1 and E[P{t)] < 2 is presented in Figure |7] 

VII. Duality relationships between stability 

REGIONS OF the MAC AND BC NETWORKS 

As discussed in the introduction of this paper, the duality 
relationship between the information theoretic capacity regions 
of MAC and BC channels in [6] provides a motivation to ex- 
plore similar results with our system model. In the discussion 
that follows, a Gaussian MAC network is defined to be the 
dual of a BC network with identical number of users if 

1) The power constraints of the MAC and BC networks 
have the same form i.e either both have average power 
constraints or peak power constraints 

2) The maximum transmit power (average or peak, as the 
case may be) of the broadcast network is equal to the 
sum of the maximum transmit powers of all users of the 
MAC network (For example, the dual MAC network 
of a two user broadcast network with power constraints 
P{t) < P has power constraints Pi{t) < Pi and P2 < 
P2 such that Pi+P2= P) 

3) Their channel states have identical statistics and the 
noise variance at all the receivers in the BC network 
are identical to the noise variance at the receiver in the 
MAC. 



4) The set of rates of the codebooks available for all the 
MAC users are identical to each other and are identical 
to the set of rates at the transmitter in the broadcast 
channel. 

We reproduce MAC-BC duality results from [6] for com- 
pleteness. 

Theorem 7.1: [Jindal et. al.] For a fixed channel state, the 
information theoretic capacity region of a Gaussian Broadcast 
channel with power constraint P is equal to the union of 
the capacity regions of the dual multiple access channel 

K 

with power constraints (Pi,P2, ■■Pk) such that "^^Pj — P- 

Furthermore, every point on boundary of the capacity region 
of the broadcast channel is the corner point of a pentagon 
representing the capacity region of some dual MAC channel. 

The fact that the boundary of the capacity region of the BC 
channel is achieved by the corner points of the MAC capacity 
regions is important in the duality between stability regions 
(below) since in our system model, only the corner points 
of the pentagon representing the MAC capacity regions are 
achievable. 

The relationship between the stability regions of MAC and 
BC networks are characterized in Theorems 17 .21 and 17.31 below. 

Theorem 7.2: The union of stability regions of all MAC 
networks that are dual to a particular BC network with peak 
power constraint is strictly contained in the stability region of 
that dual BC network. 

Proof: Refer Appendix HI] for an explanation based on 
Theorem 17.11 ■ 
However, a duality relationship in the form of equality holds 
between stability regions of peak power constrained MAC 
and BC networks of a special type - MAC networks with 
centralized power i.e the power can be shared between the 
two transmitters of the MAC network and the power constraint 
is expressed as Pi (t) + P2 (t) < P. At a given channel state 
{hi, /i2), the physical layer capacity duality result implies that, 
for every rate vector transmitted by the broadcast channel, 
there is at least one multiple access channel whose sum 
power is equal to P that can achieve that rate. Therefore at 
time slot t, powers Pi{t) and P2{t) in the centralized-power 
MAC network can be picked so as to emulate the appropriate 
'distributed power' multiple access channel to achieve the 
rate vector in the corresponding dual broadcast channel. In 
other words, the MAC channel with centralized power can 
switch between different dual MAC channels depending on 
the channel state, to achieve the rate vector transmitted in the 
broadcast network. 

For MAC and BC networks with average power constraints, 
a duality result similar to the one presented in [6] holds. 

Theorem 7.3: The stability region of a broadcast network 
with an average transmit power constraint is equal to the union 
of the stability regions of all the its dual MAC networks. 

Interestingly, the duality result holds . Proofs for Theorems 
17.21 and 17.31 are provided in the Appendices HI] and HII] . 
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Fig. 6. Stability regions of the BC and several dual MAC networks with 
peak power constraints 




Fig. 7. Stability regions of the BC and several dual MAC networks with 
average power constraints 

The plots of stability regions of BC and dual MAC networks 
with peak and average power constraints (Figures |6] and ^ 
confirm the results established in Theorem 17.21 and 17.31 

From the plot, it can be observed that in the peak power 
constraint case, the convex hull of all the MAC pentagons, 
representing time division multiplexing between the different 
multiple access channels, is a subset of the broadcast pentagon. 
Note that the convex hull of the MAC stability regions is 
different from the stability region of the MAC networks with 
centralized power; the former represents the stability region 
achieved by choosing a MAC channel among the various dual 
MAC channels, randomly. Whereas the latter stability region 
that can be achieved if a scheduler can choose a dual MAC 
channel, based on the channel conditions. 

VIII. Optimal Codebooks 

There is an interesting extension associated with the stability 
region problem in multi-rate MAC and BC networks. Consider 
a general N user MAC network. Given the number of rates 
allowed at each user i.e given |7?.i|, |7?.2|...|7^Ar|, how do 
we choose rates TZi,i — l,...Af so as to maximize the 




Number of Users (N) 



Fig. 8. Optimal codebook rate Rq{N) in a fixed codebook MAC network 
versus number of users N 



s|R;(N),») 
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Fig. 9. Optimal value of sum-rate s{Rq{N), N) in a MAC network versus 
number of users A'^ 

sum-rate point in the stability region ? The problem clearly 
has important applications in the design of communication 
systems. 

To simplify the problem, consider a MAC network with 
identical peak power constraints at all the users. Also, assume 
that all users have a codebook of rate Rq i.e \TZi\ — l,i = 
1, 2...N. Given N and Rq, the stability region of the network 
can be found using techniques presented in section [V] Let 
s{Ro,N) indicate the maximum stable sum-rate i.e 

N 

s(i?o, N) = max Rj 

i=l 

under the constraint that ...Rn) lies in the stability region. 
The optimal codebook design problem can be formulated as 
follows. 

i?2(iV) = argmaxs(i?o,A^) 
Ro 

Plots of Rq{N) (found by simulations iterating over dif- 
ferent Rq) and the coiTesponding sum-rate s{Rq{N), N) are 
provided in Figures |9] and |8] respectively. Rayleigh fading 
of unit variance, and unit power at all the transmitters are 
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Fig. 10. Average rate per user versus the number of users with the optimal 
codebook in a MAC network 

used in the simulation. Note that Figure |9] plots, versus N, 
the highest possible stable sum-rate in fixed codebook N-user 
MAC network. The increase in s{Rq{N), N) can be trivially 
explained - a sum rate of s{Rq,N — 1) can be achieved in 
a N user system, using a scheduler that simply ignores the 
A^th user The behavior of Rq{N) with increasing N is not 
known for a general channel distribution, though the plot in 
m indicates that it decreases with increasing N in Rayleigh 
fading channels. Also, the plot in Figure [TO] suggests that the 
rate seen by each user decreases with increasing A^ (This can 
also be observed simply by the fact that the rate of increase 
of s{Rq{N), N) in Figure |9] is decreasing). This suggests a 
trade-off between the the throughput of the system and the 
average rate seen by each user 

IX. Conclusions and Future Work 

A framework combining both information theoretic and 
networking points of view of MAC and BC networks has been 
presented and their stability regions have been characterized. 
The assumptions on encoding and decoding in the system 
model implied that the information theoretic capacity region 
did not present the complete picture and this led to some inter- 
esting observations. For example, we found that for networks 
with peak power constraint, the optimal scheduling policy may 
transmit at a power less than the maximum possible power in 
peak-power constrained MAC networks. 

Stability of MAC and BC networks with both average and 
peak power constraints were found. The stability regions of 
MAC and BC networks with average power constraints satisfy 
the fundamental duality property that relates the information 
theoretic capacity regions of these networks. Interestingly, the 
duality property holds in spite of restricting users to finite rate- 
sets and using a suboptimal decoding strategy. In the peak 
power constraint case, the union of MAC networks is found 
to be a subset of the corresponding dual BC networks. 

This work opens up some interesting areas of future re- 
search. Characterization of stability regions of more compli- 
cated networks such as the interference, X and Z channels, 
and relay networks can potentially provide interesting insight. 



The information theoretic capacity regions of the networks 
mentioned are not yet known; however, decoding scheme can 
be restricted successive interference cancellation and their 
stabiUty regions can be determined. Thus stability regions 
could be useful metrics for these communication networks. 
Their stability regions could possibly share some fundamental 
properties with their capacity regions like the duality property 
in MAC and BC networks. Therefore, a study of their stability 
regions may potentially provide useful information about the 
nature of their capacity regions. Another area of future work 
is the design of optimal codebooks for these networks dis- 
cussed in Section rVIIII While Section [Villi discussed the fixed 
codebook version of the problem, an interesting issue is the 
loss in performance by restricting users to a finite number of 
codebooks in comparison to the information theoretic capacity 
region. In other words, a characterization of the number 
codebooks that should be used so that the sum-rate achieved is 
'close' to the maximum sum-rate in the information theoretic 
capacity region is an important related problem in the design 
of communication networks. 
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Appendix I 
Proof of Lemma 1 

Proof: The fact that a channel state can be uniquely 
mapped to a single rate in the optimal policy in networks 
with average power constraints can be proved for general N- 
user MAC and BC networks with multiple rate-sets but such 
a proof is complicated. The idea is conveyed in the proof of 
the result for a fixed codebook point-to-point link. 

Consider a point-to-point link with average power con- 
strained by E[P{t)] < P. The link has a fixed codebook 
of rate Rq. The channel is Gaussian with gain h{t) which 
takes values from a continuous state space. At a given state, 
h, the only parameter that affects the average rate and power 
is the probability of transmission at this state p{h). Note that 
the effect of any algorithm, in terms of average rate and 
power consumed can be achieved by an equivalent randomized 
algorithm with the appropriate 'weight' function p{h). The 
problem of finding the stability region of this network is 
simply equivalent to 

max(/i) / RQp{h)f{h)dh 
P J 

s.t J P{h)p{h)f{h)dh < P 

where P{h) is the power required to transmit a codeword at 
state h. Note that P{h) is monotonically decreasing in h. 

We intend to prove that in the optimal policy a given 
channel state is mapped to a single rate vector in {i?o,0} 
or equivalently, that p{h) G {0, l},V/i 

Let the optimal rate allocation policy be p*{h). Contrary to 
our claim, let p*{h) take values in (0, 1) for certain G R. 
Now, consider an interval (x, j/) along the positive axis on the 
real line such that p{h) ^ {0, 1}, V/i G (x, y) and /J f{h) > 
. Note that if we cannot find such an interval, then setting all 
values of p'{h) = [p(/i)J is also optimal and satisfies the 
given property. We now find an alternate scheduling policy 

which performs better than p*{h) as a contradiction. 
We know that since p*(h) < 1, V/i G {x, y) 



p*{h)P{h)f{h)dh < / P{h)f{h)dh 



3ze{x,y)s.t r p*{h)P{h)f{h)^ r P{h)f{h)dh 

P{h)p*{h)f(h)dh = J' - P*{h))f{h)dh 

Since P{h) is decreasing in h, we have 

p{z) ( p*{h)f{h)dh< r p{h){i-p*{h))f{h)dh 

J X J z 

and 

' P{h){l-p*{h))f{h)dh < P{z) [\l -p*{h))f{h)dh 



Combining the above two inequalities, we have 



p*{h)J{h)dh< / {l-p*{h))f{h)dh 



p*{h)f{h)dh< / f{h)dh 



(5) 



Now, we construct a new scheduling policy as follows 

r p*ih) ifhi{x,y) 
pi{h) = < if he[x,z) 
[ 1 ifhG [z,y] 

Equation |4] implies that the power consumed with pi{h) is the 
same as the power consumed with p*{h). Equation |5] implies 
that the average rate achieved with pi{h) is greater than that 
achieved with p*{h), contradicting the assumption that p*{h) 
is optimal. The result is hence proved. ■ 

Appendix II 
Proof of Theorem |7.2| 

Proof: We observe the channel state space {hi{t), h2{t)) 
has identical statistics for the broadcast and the dual MAC 
channel. Therefore, if we show that, for a given channel state, 
all rate vectors that can be supported by the MAC channel 
can be supported by the dual broadcast channel, then the 
stability region of the multiple access network is a subset of 
the network capacity region of the broadcast network (since 
the BC network can follow use the appropriate scheduling 
strategy and achieve the corresponding rate vector) . This is 
precisely the result stated by the physical layer duality result 
(Theorem IVIII i. and hence the result follows. 

In fact, the union of MAC stability regions is a proper 
subset of the BC stability region. To see this, we perform the 
following thought experiment. Considering a two-user fixed 
codebook broadcast channel with a peak power constraint 
P whose channel gains {hi^h^) can, take two states. Si ~ 
(/iii,/i2i) and S2 — (/ii2,^22) with equal probability (of i) 
such that 

hi2 = /121 < hii = /122 



and 



log(l + hl^P) = i?o 



where Rq is the rate of the codebook in the BC channel. 
Therefore, we have : 

log(l + hl^P) < Ro 




BC stability region 



(-R.0,0) 

Union of dual MAC stability regions 



Fig. 11. Stability regions of BC and dual MAC networks with 'ON-OFF' 
channels 



In Other words, there two possible channel states. In each 
channel state, one of the two channels is ON and the other is 
OFF. The corresponding set of supported rates are {(i?o,0)} 
and {(0, Rq)} (assuming additive white Gaussian noise of unit 
variance at each receiver). Using the result of equation [T] the 
stability region of this BC network is a square in the Cartesian 
plane (see Figure fTTTl formed by the points (0,0), (i?o,0), 

(0, Ro), (Ro, Ro)- 

Among the set of all dual multiple access networks with 
power constraint Pi and P2 satisfying P1+P2 — P, transmis- 
sion is possible only if Pi = P or P2 = P. In other words, if 
Pi < P and P2 < P, no transmission is possible. In the MAC 
network with power constraints (P, 0), the sets of supported 
rates are {{Rq, 0)} for channel state ON and (0, 0) for channel 
state OFF. The stability region is therefore a line joining 
(0,0) and (-|i,0). The union of stability regions of dual MAC 
channels is thus, the union of two line segments along the 
axes forming the adjacent sides of the square representing the 
BC stability region. Clearly, this is a proper subset of the BC 
stability region. This proper subset relationship is confirmed 
by results in Figure |6] 



We construct the scheduling policy as follows. Since fl 
lies in the capacity region of the broadcast network (for 
h), we know from Theorem 17.11 (Section IVIII ) that there 
exists a power allocation {Pi,P2) in the MAC network we 
have considered, such that Pi + P2 = P{h) and fl{h) can 
be transmitted successfully. Note the capacity region of the 
broadcast channel is achieved by the corner points of the 
pentagons representing physical layer capacity regions of the 
dual MAC channels. Therefore, in spite of disallowing rate- 
splitting, we can find {Pi{h), P2{h)) in the MAC channel so 
that fl{h) can be transmitted. Also, for a given h. Pi and 
P2 are purely functions of the rate vector jl{h) transmitted 
by the broadcast channel. Let Pi{h) — f{h,fl{h)) and 
P2{h) = g{h, j2{h)). Now the scheduler in the MAC networks 
transmits fl{h) with the two users using powers f{h, jl{h))) 
and g{h,jl{h)). Clearly this MAC network can stabiUze any 
arrival rate stabilized by the broadcast network. Also, since 
E[Pi{h)] + E[P2{h)] = E[Pi{h) + P2ih)] = E[P{h)], and 
therefore E{g{h, il{h))) + E{f{h, fl{h))) = P, the power 
constraints of the MAC network satisfy the duality constraint. 
Therefore, the broadcast network stability region is a subset 
of the union of stability regions dual multiple access channels 
with average power constraints. An analogous argument can be 
used to show that the union of stability regions of dual MAC 
networks is a subset of the stability region of the broadcast 
network. Thus the stability region of a broadcast network with 
average power constraint is equal to the union of stability 
regions of dual MAC networks. ■ 



Appendix III 
Proof for Theorem I7.3I 

Proof: In a BC network with average power constraints, 
the optimal policy schedules a rate vector /l(/i) when the 
channel state is h. Let the power transmitted at channel state h 
be P{h). We provide a MAC network and a scheduling policy 
for this network that achieves the following: 

1) The average rate vector as a result of the scheduling 
policy is equal to the average rate vector of the BC 
network 

2) The average power expended by the users in the MAC 
network satisfy the duality property i.e, the sum of the 
average powers consumed is equal to the average power 
expended in the BC network. 

Existence of such as scheduling policy ensures that the stabil- 
ity region of the broadcast network is a subset of the union of 
stability regions of dual MAC networks. 



